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ABSTRACT 

Scalar field theory at finite temperature is investigated via an improved renormal- 
ization group prescription which provides an effective resummation over all possible non- 
overlapping higher loop graphs. Explicit analyses for the A(/)^ theory are performed in 
d = 4 Euclidean space for both low and high temperature limits. We generate a set of cou- 
pled equations for the mass parameter and the coupling constant from the renormalization 
group flow equation. Dimensional reduction and symmetry restoration are also explored 
with our improved approach. 
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I. INTRODUCTION 



In recent years, intensive efforts liave been devoted to quantum fieid tlieory at finite 
temperature, a subject witli wide applications in areas sucli as tlie evolution of the early 
universe and its cosmological consequences [1] , the deconfinement phenomena and forma- 
tion of the quark-gluon plasma [2] , and the critical behavior of condensed matter systems 
near the phase transition. To investigate these issues, one generally utilizes the finite- 
temperature effective potential approach in the spirit of the perturbative loop expansion. 
However, in the high temperature limit, perturbation theory becomes unreliable [3] since 
the presence of infrared (IR) divergences may destroy the correspondence between the 
expansions of loops and the coupling constant. Certain higher loop contributions such as 
the "daisy" and "superdaisy" diagrams that contribute to the same order in the coupling 
constant must also be incorporated [4] for computing the critical transition temperature 
Tc and determining the nature of the phase transition. In gauge theories, it has been shown 
that the "hot thermal loops" need to be resummed in order to obtain a gauge independent 
gluon damping rate [5]. 

However, since resumming the multi-loop contributions to arbitrary orders proves to 
be nontrivial, various methods have been proposed to carry out the resummation: In [4], 
the use of a gap equation was first discussed; in [6] a renormalization group (RG) with 
T, the temperature, as the fiow parameter was used; and in [7] a self-consistent Hartree- 
Fock formalism is presented. An "environmentally friendly" renormalization prescription 
for interpolating effective finite-temperature theories in different regimes has also been 
proposed [8]. Although all of these techniques yield the same finite-temperature effective 
propagators in the leading order, they differ in the subleading correction which can affect 
the nature of the phase transition. 

In this paper, the methodology adopted for investigating the scalar field theory at finite 
temperature is based on the use of RG constructed from the Wilson-Kadanoff [9] blocking 
transformation in the Euclidean formalism. Unlike [6] , our RG arises from the arbitrariness 
of the internal blocking scale k, and not the external temperature parameter T. The 
formulation not only takes into consideration the dominant higher loop diagrams without 
resorting to the complicated analytical order by order resummation, it also characterizes 
the flow pattern of the theory for arbitrary T as well as the momentum scale k which 
is chosen to be zero in the methods described above. Our approach is analogous to the 
series of works by Tetradis and Wetterich [10]. However, instead of using a smooth 
momentum smearing function which leads to an integro-differential RG equation, our sharp 
momentum cut-off [11] yields a RG flow equation which takes on the form of non-linear 
partial differential equation. Besides the advantage in performing numerical computation, 
it also offers a more lucid physical interpretations. 

Consider the following scalar lagrangian: 



At finite temperature, the RG improved blocked potential Up^k{^) associated with the 
blocked field ^{x) is characterized by the following differential RG equation [11] : 




-(9^0)2 + ^(0). 



(1.1) 




(1.2) 
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The above equation is obtained by first computing the finite-temperature blocked potential 
to the one-loop order followed by a RG improvement to take into account all possible non- 
overlapping daisy, superdaisy and higher-loop diagrams which strongly modify theory in 
the high T and small k limits. We also see that in this RG construction the contributions 
of a particular mode which has been integrated out are naturally retained for the integra- 
tion of the next, and therefore the interactions among the modes are properly taken into 
account. In addition, during the course of mode elimination, irrelevant operators defined 
with respect to the ultraviolet (UV) fixed point continue to be generated and their ef- 
fects incorporated as well by (1.2) [12]. Furthermore, our flow equation also offers insights 
into the existence of high temperature dimensional reduction as well as the restoration of 
symmetry accompanied by the disappearance of imaginary contribution to C//3,fc($), as we 
shall see later. To describe the full theory, however, one also needs to consider the effects 
of wavefunction renormalization constant Zp^i.{^) as well as the higher-order derivative 
terms in the blocked lagrangian. Nevertheless, the use of (1.2) is justified in the IR limit 
of the four-dimensional theories since the other contributions are only of higher order. 

The organization of the paper is as follows: In Sec. II the formalism of the blocking 
transformation at finite temperature is briefiy reviewed using the scalar Xcf)'^ theory as 
illustration. The finite temperature RG flow equation (1.2) with underlying 0(3) symmetry 
is constructed and compared with the simple one-loop independent mode approximation 
(IMA) . In Sec. Ill we concentrate on the low T regime where thermal effects are negligible 
and compare the 0(3) thermal blocked potential C//3,fc($) with the 0(4) symmetric Uk{^) 
at T = 0. The high T limit of the theory is investigated in Section IV. The criterion for 
dimensional reduction is explicitly deduced from (1.2). In addition, the behaviors of the 
scale-dependent thermal mass and coupling constant ^ and Xp^u, are studied with a set 
of coupled equations. It is shown that Xp^k decreases with increasing T and approaches 
a constant in this regime. The complementary relationship between the internal scale 
k and the external parameter T is also demonstrated. The phenomenon of spontaneous 
symmetry breaking at T = and symmetry restoration above Tc is discussed in Sec. VI. 
The phase boundary is determined from (1.2) by requiring the imaginary part of Up^ki^) 
to vanish at Tg. The value obtained in this manner is compared with that obtained in the 
one-loop approximation without the imaginary sector. Sec. VI is reserved for summary 
and discussions. We collect in the Appendix the details of extracting the leading order 
thermal contributions from the integrals encountered in Sec. IV. 

II. FINITE-TEMPERATURE RENORMALIZATION GROUP 

In this section, we briefiy review the finite temperature RG formalism used in investi- 
gating the thermal behavior of the theory. The concept of the RG is based on the notion 
that in certain physical processes only a particular range of modes in the momentum de- 
composition of the field (f){x) will be relevant, and it is often desirable to eliminate the 
irrelevant modes to which the physics is insensitive. The reduction of degrees of freedom 
is then compensated by a readjustment of the parameters in the coupling constant space 
of the lagrangian. Therefore, instead of using the original field variable 
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where /3 ^ = T and ujn = 2nn/P denotes the Matsubara frequency, we define the coarse- 
grained blocked field as [11]: 



$(x) =Tj^^ dy^ J d^ypfc(x - y^iy) 



(2.2) 



via an 0(3) invariant smearing function pfc(x). Since the low energy physics is unaffected 
by the modes above the "blocking scale" k, we shall for simplicity choose 



Pk 



(x)= / 

^|p|<A; 



(27r)= 



(2.3) 



or pk{p) = 0(fc — IpD- In other words, k acts as an upper cut-off for the modes which 
are to be retained. Notice that in this formulation, the S- function associated with the 
imaginary-time variable r is given by 

oo 



')4 ^ 

' 'n — — rv 



(2.4) 



Since the wavefunction renormalization constant Zk{^) yields only a minute correction 
to the anomalous dimension, we shall for simplicity set Zk{^) to be unity. In addition, 
by neglecting the higher order field derivative terms, one may simply choose the static 
limit $(x) = $ = const, for computing [/^ fc($). Upon a simple Gaussian integration, the 

one-loop contribution to the finite-temperature blocked potential C//3^/e($) becomes 



' n — — r 



In 




(2.5) 



P^/p'^ + V"{^) + 21n 
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where the prime on the p integral indicates an integration subject to the constraint k < 
p < A, with A chosen as the UV regulator for the three-momentum integration. For 
F($) = //^$^/2 -|- A$^/4!, by demanding the renormalized parameters to satisfy 



a$2 



4>=/3- 



=fc=0 



=fc=0, 



(2.6) 



the resulting Uj3^k{^) takes on the form: 
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1/2 
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,~2 ^2^^ ik + ^k^ + iii + XR^y2- 

+ [Pr + ^>^R^ ) In 



(2.7) 
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27r2/? 



which in the hmits of vanishing and fc, reproduces the usual effective potential [13]. 
Therefore, one sees that consistency with the 0{4) invariant theory in the T = limit 
requires a complete summation over the Matsubara frequencies ujn = 2nn/P for all n, 
i.e., Pfc(p) must be independent of Un- It is also evident from (2.6) that temperature- 
independent subtractions alone are suflBcient to remove the cut-off dependence. 
Differentiating (2.5) with respect to the arbitrary scale k leads to 



= -^Vk' + V"{^)- T^lnW - e-^v/^^^^l . (2.8) 
ok 47r^ 27r^ L J 

This RG equation is obtained in a manner in which each individual mode is integrated out 
independently by neglecting the systematic feedbacks from the high modes to the lower 
ones during the elimination. The parameters for this independent-mode approximation 
(IMA) scheme are denoted by a tilde. 

On the other hand, instead of integrating out all the modes from A to all at once, 
one may first divide the integration volume into a large number of thin shells of small 
thickness Ak. By lowering the cut-off infinitcsimally from A — A — Ak until A = /c is 
reached, we arrive at the following RG equation governing the flow pattern of the improved 
finite-temperature blocked potential Up^ki^)'- 



j^dUp^__k^ I ^ , _T— In 



This non-linear partial differential equation establishes a smooth connection between the 
small- and large-distance physics at finite temperature. Moreover, (2.9) systematically 
incorporates the contribution of a particular mode for the elimination of the next. In 
fact, (2.9) takes into account all possible non-overlapping daisy, superdaisy and higher- 
loop diagrams which significantly modify the theory at high T and small k limits. The 
summation over the multi-loop graphs from the basic one-loop structure is depicted in 
Fig. 1. In other words, what we have here is a "physical" one- loop diagram characterized 
by the dressed vertices. One may certainly improve (2.9) by including higher "physical" 
loop contributions. However, we argue that since each loop integration is multiplied by 
a factor k = Ak/k, a term of order m in loops will be suppressed by k"^ in the limit 
Ak 0. Therefore, the use of (2.9) is justified. What we have neglected in (2.9) are 
the overlapping higher loop diagrams. However, as shown in [4], if one considers the A^- 
component theory and takes the limit — > oo, these overlapping graphs will be suppressed 
by an extra factor N compared with the corresponding non-overlapping contributions at 
the same loop order. The effects of these overlapping graphs for the one-component case 
is currently being explored. 

Eq. (2.9) implies the following RG fiow equation for the effective scale-dependent 

thermal parameters = Up^^{Q), Xp^k = t^^^fe(O) and gp^k = uj^%{0)- 



dk 87r2 



_|_ ^2 ^^ 



coth — !^ , (2.10) 
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and 
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dX 



I3,k 
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167r2 k^ + ^2 , 
^9f3,k k^ 



87r2 




2 

/3,fe 



T{e^V''^+^0'k _ ly 



(2.11) 

Notice the contribution to the flow of X/^^k from the sixth-order coupling constant g^^k- The 
presence of such higher order corrections is the natural consequence of blocking transforma- 
tion through which the number of effective modes is reduced at the expense of generating 
a more compHcated effective action. 

While Up^k{^) associated with the IMA is solved analytically in (2.7) apart from the 
integral containing thermal effects, the RG improved Up^ki^) can only be obtained by 
solving (2.9) numerically. In Figs. 2 and 3, the evolution of Up^ki^) and Up^ki^) foi" 
T = and T 7^ are traced. We notice that the two potentials agree reasonably well 
for k near the cut-off A and begin to deviate as k is lowered. In the deep IR limit where 
A; = 0, appreciable difference between t/^_fc($) and C/^,fe($) is seen. Such a difference can 
be understood by noting that 



k^ 



47r2 



27r2 



In 



1 _ e-V='+^^U(*) 

1 _ g-/3Vfc2+V"(*) 
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exp{-(3{^k^ + V"{<5)- Jk^ + U'^,ki^))} 



27r2 
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1 - exp{-f3y/k^ + V"{^)} 



2T 
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X 
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4 V P + V"{^) 



+ 



(2.12) 

Thus we see that higher loop corrections will continue to pile up as k is lowered or T is 
raised. By comparing Fig. 2 and Fig. 3, we see that the inclusion of the thermal effects 
improves the agreement between the RG and the IMA results. In Sec. IV we shall see that 
the high T behavior of the theory is strongly modified by the higher loop contributions. 

In Fig. 4 the evolutions of ^ and A/3,fc with k for three different values of T are 
depicted. While ^ increases with T, Xp^k is shown to decrease with rising T. On the 



6 



other hand, when the role of k is considered, its influences on /i^ ^ and A^^^ are completely 
opposite to that of T. One therefore observes an interesting competition between the exter- 
nal parameter T and the internal momentum scale parameter k from which RG is defined. 
By external parameter we mean a "physical" quantity which has direct physical impact on 
the system and can be measured experimentally. On the other hand, an internal scale is 
a "fictitious" scale which we choose to characterize the theory. In this case, k is chosen to 
provide a separation between the high and the low modes. However, the precise definition 
of "high" and "low" modes is dependent on the energy regime under investigation. A more 
complete discussion concerning the behaviors of the thermal parameters for different k will 
appear in the later sections. 



III. LOW TEMPERATURE LIMIT 



We first consider the low temperature limit T « ^k"^ + C/^'^ where, upon dropping 
the subscript j3, (2.9) is reduced to 



" dk - 47rW^ 



k 



(3.1) 



This "zero" temperature limit can be compared with the RG flow equation associated with 
t/^($), the blocked potential for d = 4 derived in an 0(4) symmetric manner [14]: 



. . k^ 

kdtUkm = -^ln( 



(3.2) 



Note that quantities derived from the 0(4) symmetric t^^($) shall be distinguished with 
a caret, and double carets for the corresponding IMA scheme. 

Since the underlying symmetry of the finite temperature RG equation (2.9) is 0(3), 
it is instructive to compare its low T flow patterns with (3.2) for both RG improved and 
IMA prescriptions. Eqs. (3.1), (3.2) along with (2.8) lead to the following RG coefficient 
functions for the running mass parameter and the running coupling constants: 



(i2 = k 



P2 = k 



(3, = k 



132 





Xr k^ 


dk 


Sn^ ^k^ + 




Xk k^ 


dk 






Xr k^ 


dk 




m 




dk 


k 



[0(3), IMA] 
[0(3), RG] 
[0(4), IMA] 
[0(4), RG] 



(3.3) 
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p4 = k 



(3^ = k 



(3^ = k 



/?4 = k- 



dXk 
dk 

dXk 
dk 

dk 



3A| 



k^ 



3A| k^ 
k^ 



9k 



3Ai 
167r2 



3A| 



k^ 



9k 



k' 



dk 



{k'^ + plf 



167r2 fc2 + ^2 



[0(3), IMA] 
[0(3), RG] 
[0(4), IMA] 

[0(4), RG]. 



(3.4) 



Notice that in the RG improved 0(3) and 0(4) schemes for A^, the influences from 
is incorporated. For k"^ » and k^ » fi^, one finds agreement in the leading order 
behavior of (34 and P4 for the running couphng constant, as can be seen from (3.4). How- 
ever, for the running mass parameter, and P2 differ by a factor of two in this Hmit. 
Such a discrepancy can easily be understood by noting the difference in the symmetries of 
the underlying RG constructs, namely, 0(3) for the former and 0(4) for the latter. The 
manner in which the original bare theory is approached can be elucidated by examining 
the simple one-loop structure of the bare mass parameter: 



2 ~2 



A 



R 



327r2 

Xr 
327r2 



2A2-/i|,ln(^)+/i|(l-ln2) 

I^R 

H'R 



:0(3) 



0(4), 



(3.5) 



where A and A are, respectively, the three- and four-momentum ultraviolet cut-offs. There- 
fore, in the low-temperature limit with k'^ » Up^ki^) can be transformed to Uj^{^) 
with the approximate scaling relation k k/\^. The reason for the existence of scaling 
can be explained by noting that the one-loop zero temperature contribution to Up^ki^) 
from (2.5) can be rewritten at /c = as: 



where we have used 



(3.7) 



which holds up to an ^^-independent constant. With the coordinate transformation 

Po = P sin^ 
p = p cos9 



(3.8) 



8 



and imposing the same cut-off regulator A for po such that < < 2A^, the above 
expression becomes: 



f7r/2 

47r2 



1 



A 



dpo In 



COS^^ 



1 



167r2 

_ 1 
~ 2 



V2A 



dp In 

d'^p 



In 



dp p^ln [p2 + V^"($) 



p' + V"{^) 



(3.9) 



Therefore, wc sec that in the UV regime where A; ~ A, with the rescahng A = A/y/2 — > oo, 
the results derived from 0(3) and 0(4) RG schemes agree with each other. 

On the other hand, in the IR limit where « /j,^ and k^ « ft^, we have 



P2 = - 






[0(3), IMA] 






Xkk^ 


1 /c^x 


[0(3), RG] 


k = - 




1- +••• 


[0(4), IMA] 






k^^ 1 




[0(4), RG], 


167r2/i2 I 



(3.10) 



and 



/34 



167r2/i| 



167r2/i|, 

3A|fc4 
fc /I 



2/i2 

3P 



Ay 

2fc2 



+ 



+ 



+ 



167r2/i2 



+ 



[0(3), IMA] 
[0(3), RG] 
[0(3), IMA] 
[0(4), RG]. 



(3.11) 



Comparing (3.10) and (3.11), one sees that the running parameters increase more rapidly 
for the 0(3) symmetric case near A; by a factor of //fe//c in the IR regime, and the 
scaling relation between k and k for transforming t/g^/-($) and ^7^($) in the UV regime 
is clearly lost. This is due to the fact that (3.7) no longer holds for finite k and that a 
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transformation similar to that in (3.9) for finite k ceases to exist. The numerical results 
for the flow equations in (3.3) and (3.4) are illustrated in Figs. 5 and 6. For the mass 
parameters shown in Fig. 5, the flow begins at A = A/\/2 where all schemes take on the 
same negative bare mass parameter /x^. We have chosen /x^ < although in principle it 
can also be positive. However, the physics we are interested lies in the IR regime which 
is very far from the cut-off. Therefore, the low energy physics is influenced by the sign 
of the renormalized mass parameter but not of /U^. The scaling relation between the 
0(3) and the 0(4) schemes in the large k limit is seen by noting that the two curves are 
approximately parallel to each other initially. One also observes that both RG improved 
prescriptions yield the same mass parameter at = 0, independent of the underlying 
symmetry of the RG. Such a behavior is also observed within the IMA schemes. This 
symmetry-independent behavior is again due to the transformation (3.9) which leads to 
identical one-loop effective potential. The same argument can also be made for the RG 
improved prescriptions by simply replacing ]/"($) in (3.9) by the corresponding ,^{^). 
This also explains the necessity of having a more rapid convergence in the deep IR limit 
in order that both RG schemes yield the same improved parameters. 

Another notable feature arises from comparing IMA with the RG approach, where 
we see that by choosing an initial condition /i^ < 0, /x|^o/A*fi ~ 4.6 for A = 7. Such 
an enhancement in the mass parameter originates from the systematic accumulation of 
higher loop contributions in the course of repeated blocking transformation which brings 
the theory to the IR regime. On the other hand, if //^ > is chosen, the RG improved 
/i^^o ^^^^ ^® comparable to the IMA result nj^. 

The remarks on the running of the mass parameter also apply qualitatively to the 
flow pattern of the coupling strength which is displayed in Fig. 5. However, contrary 
to the former, there exists only a minute difference between the RG and IMA results. 
Therefore, the coupling constants obtained with various approaches near zero temperature 
are approximately equal. As we shall see in the next section, this is not the case for the 
large T limit, a regime where the IMA scheme becomes unreliable. 

IV. HIGH TEMPERATURE LIMIT 



We consider next the high temperature limit where T >> + t^^'jf Since in this 

temperature range the second term in (2.9) dominates, if one neglects the first term in 
(2.9) entirely, the RG equation is simplified to 



^ k^ rfc^+Cg^ 

dk 47r2 Lp + t/^/^(o) 

or 

dk 47r2 Lp + ^^(o) 



(4.1) 



(4.2) 



where Uk{^) — (^Up^ki,^) and the dots denote differentiation with respect to i> = v^*^? the 
new field variable defined in d = 3. One therefore concludes that the theory in this limit has 
undergone a dimensional reduction (DR) with temperature being completely decoupled. 
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leaving an effective three-dimensional theory described by (4.2). This is precisely the flow 
equation one would obtain for a three dimensional theory at T = 0. In other words, the 
high T behavior of the theory in d = 4 corresponds to that of (i = 3 at zero temperature. 
Physically this phenomenon can be explained by noting that the high temperature limit is 
characterized by the shrinking of the "imaginary time" dimension having a period (3. This 
in turn implies a suppression of the n ^ non-static modes in the Matsubara summation, 
thereby giving vanishing contribution to C/(3^fe($). By neglecting the non-static modes 
completely, the remaining static sector is what the three-dimensional theory parameterized 
by Uk{^) describes. Notice that the coupling constant would become Xn = XjiT carrying 
the dimension of mass. 

The above analysis demonstrates that such a dimensional reduction takes place only if 
(fc^ -I- Uj^ f^{^))/T'^ 0, where C/^'^(0) is the thermal mass parameter /i| j^. The condition 
also implies a small momentum scale k « T. However, as noted in [15], DR strictly 
does not takes place in the infinite temperature limit because fXp acquires a T^-dependent 
correction in the leading order which renders the ratio U'^ /T^ finite for all T. In Fig. 7, 
we compare the parameters obtained from (4.2) with that of the full RG equation (2.9). In 
both mass and coupling constant, we see that the two results differ by a finite constant gap 
which persists to arbitrary large value of T. This indeeds supports the conclusion drawn 
in [15]. Nevertheless, for sufficiently large T, the theory exhibits a "partial" dimensional 
reduction since the minute difference can be neglected. We also notice that the value of 
at T = predicted by the DR prescription corresponds to the bare coupling constant 
Xb > Xr. One must remember, however, that the results obtained from the dimensioanlly 
reduced prescription for low T are not reliable. 

Another interesting issue one may explore in the high temperature limit is to examine 
the behavior of the effective thermal coupling constant Xp^k [16] and demonstrate the 
ineffectiveness of the one-loop IMA scheme. To see how perturbation theory breaks down, 
we notice first that from the IMA, the effective scale-dependent thermal parameters can 
be written as [11]: 




(4.3) 



and 




(4.4) 



oo 



Z 




-1 



= Xk 



87r2 



dz 



a^){z'^ 6^) 
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where o? = P^{k^ + and 6^ = and the last integral expression in (4.4) is obtained 
via an integration by part. In the above, 



A 



R 



64712 ^FT/^ 



k + ^/W^^p^ 



(4.5) 



and 



2.1], / A:[2A: + (2F + /l|)(F + yj22^)-i/2] 



-H 



k + ^/l^T^- 



(4.6) 



represent the running parameters at T = 0. The integrals appearing in (4.3) and (4.4) 
from finite temperature contribution can be approximated in the limits of large or small a 
and b. The details arc provided in the Appendix following the classic treatment by Dolan 
and Jackiw [4]. With (A. 15) and (A. 34) in the Appendix, the one-loop approximation for 
the thermal parameters in the small a and b limits take on the forms 



~2 _ ~2 , 



A 



A 



R 



Rl^R 



24 
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3XrP XriIrT 
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A'2 



AttT 



47r2 



sin 



+ 7-1- tanh 



k 



(4.7) 



and 



^I3,k — ^k 



3A2jT r TT 



— — sm 



-u ^ ) 



+ 



3A^fc 
327r2T 



3A| 
167r2 



In I 



f^^?5^)+7+ tanh-f^i=) 



(4.8) 



Eq. (4.7) shows that at sufficiently high T the mass parameter grows quadratically 
with T. However, the presence of — 3Aij/c2/167r^ tends to decrease /i| ^. Again, we observe 
a competition between the internal scale k and the external parameter T. 

On the other hand, for large T, one naively expects the negative^ linear T-dependent 
term to dominate (4.8) giving rise to a vanishing or even negative Xfj^k- However, it is 
a well known fact that the correct high T behavior of A/3 fc cannot be accounted for by 
the simple one-loop result (4.8), and that when higher loop contributions such as daisy 
and superdaisy graphs are incorporated, the linear term will be suppressed. Since it is a 
nontrivial task to resum these higher loop effects, various attempts based on the use of 
gap equations [17], the Schwinger-Dyson equation [18] or RG have been made to provide 
an effective resummation. However, the details of resummation have raised some concerns 
[19]. We shall now illustrate how the RG equation obtained in (2.9) can be utilized to 
addressed these issues. 



12 



As explained in sec. II, a resummation of all possible non-overlapping diagrams can 
be achieved with our RG equation (2.9). By comparing (2.9) with (2.8), we see that the 
difference between the IMA and the RG prescriptions is due to the fact that in the latter, 
instead of using the /c-independent V{^), Up^ki^) is employed to give the scale-dependent 
thermal parameters /i| ^ and Xf^^k which are then used in the evaluation of Up^k-Aki^)- By 
iterating this procedure, higher loop contributions are automatically taken into account. 
Guided by this logic, we proceed to improve (4.7) and (4.8) by replacing the right-hand- 
side of the expressions by the effective /i| and Xp^k- This leads to the following set of 
two coupled equations in the limit of small a = l3{k'^ + ji^^ /j)^/^: 



Xfj.k.T^ 



,,2 _ ,,2 , ^/3,fc 



24 



+ 



167r2 



In 
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TT 
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and 



Xp,k = Xk 



^'7^^fJ'l3,k 



^^%k 
167r2 



In 



TT 

2 



k 



A-2 



AttT 



sin 

y + i^/3,k 

^ +-f+ tanh"^ 



+ 



327r2T 



(4.10) 



These coupled equations are actually inferred by (2.10) and (2.11). The only difference is 
that we have truncated the higher order contributions in (4.10) and (4.9). Notice that they 
are slightly different from those described by Chia [17] in that the continuous feedbacks 
from n'p ^ to Xf3^k and vice versa are systematically incorporated, hence leading to a more 
accurate determination of the high temperature behavior of the theory. 
Taking the k = limit for simplicity, we have 



f4 



A/jT^ X(3 upT Xf3 li} 



24 



87r 



167r2 



In 



(4.11) 



and 



Xf3 = Xr 



3X}T 



3X} 



IGtt//^ IGtt^ 



In 



4nT 



) +7 



(4.12) 



where the subscript k is dropped for brevity. From (4.11), we notice the quadratic T 
dependence of the thermal mass parameter /x^. To probe the high T behavior of A/j, 
however, would require a more precise determination of the ratio TZ — jip/T present in the 
equation for A^. In the one-loop approximation, the use of up = nn for (4.12) leads to a 
rapid linear decrease of A^ with T, signalling the breakdown of perturbation theory. By 
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including the dominant XpT'^ /2A term which gives as T — > oo, we obtain 

from (4.12) the relation 

A^ = A, + ^Af + ..., (4.13) 

which for = 0.1 takes on the value A/j = 0.09186. Further refinement with TZ = 
A^ - 3A/3/487r yields 

or A/3 = 0.09160. This is in accord with that illustrated in Fig. 8 for the T dependence of 
the thermal parameters up and Xp. Thus, we conclude that at very high T, Xp approaches 
a constant nonvanishing positive value. In Fig. 7, the infinitesimal decrease of Xp in the 
high T regime is due to the nature of the coupled equations in which a particular value of 
Xp is first used in (4.11) for deducing TZ which is subsequently inserted to (4.12) for the 
determination of a new improved Xp. In general, the more accurately TZ is, the better the 
agreement with that generated by (2.9). We emphasize here that contrary to the claim 
in [19], Xp can never increase with T because of the negative sign associated with the 
temperature-dependent term in the full RG flow equation (2.9). The inclusion of higher 
loop effects can only modify the leading T-dependent behavior, but not the sign. When 
RG is invoked to study A/3 jt, it is crucial to keep in mind that there are three dimensionless 
combinations b — k/T, h' — j^u/k and h" — [IrIT present in the flow equation with h' and 
h" having opposite effect to the running of Xp^^. Erroneous use of h" in the RG analysis 
can lead to the wrong claim that A/3 fc rises with T. By employing the correct choice 6', 
one will reproduces the standard result that Xp^k increases logarithmically with k. 

We also notice that the mass parameter /i^^ obtained with the RG method is slightly 
smaller than that of the one- loop result. The reason again is due to the continuous feedback 
of the higher loop effects as well as the use of improved value for the thermal coupling A^ 
which is smaller compared to Xr. 

Since the scales k and T enter in a complementary manner, in the regime where 
k » T one expects the thermal effects to be suppressed. This can be seen by noting that 



and 



/^L = Al + x^e-«[l + ^ + g + ... , (4.15) 



A/3, A: — A/e — ^ 



"W^ p-a p 1 

'^^'^ ^ W + a-^{a'-b') + --- (4.16) 



87r2 , 

with the help of (A. 19) and (A. 37). Indeed, the thermal parameters are suppressed expo- 
nentially by a factor e~". When the blocking scale k coincides with the UV cut-off A, all 
thermal contributions must vanish since for k = A the original temperature-independent 
bare theory is recovered. Equivalently, one may say that the counterterms are independent 
of temperature. This implies that we must have e~^^ — > in principle. Correspondingly, 
we work within a temperature range which lies sufficiently far from A to ensure that nearly 
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all particles will have a momentum below A according to the Bose-Einstein distribution. 
The suppression of thermal effects for large pk can indeed be seen from Fig. 4. 



V. SPONTANEOUS BREAKING AND RESTORATION OF SYMMETRY 

Turning to the scenario of spontaneous symmetry breaking (SSB) with /j,^ < 0, the 
one-loop finite temperature blocked potential becomes 



= f - i^) + if 



Xr 36jl% + 84XRfi%M^ +25X%M 



[2 nj4 



2567r2 



+ 



^ I -2k (2^ + il'R+ ^ Afl^^) (k' + fl% + ^Xr^') 



{f,l + XRMy2y 

1/2 



+ (/ifl+^Aii$2j'ln 



p2A;2 + fj,l + Afl$V2 + 2A;yP + + Xr^^/2^ 
fil + ~XRMy2 



+ 



27r2/3 J, 



where we have chosen the following off-shell renormalization conditions used in [13]: 



(5.1) 



ill = Re 



$=/3-i=fc=0 



< 



A 



R 



(5.2) 



a$4 



$=M,/3-i=fe=0, 



where M =< $ > is the nonvanishing vacuum expectation at which Up^ki^) is 
Since /i^ < 0, it can no longer be interpreted as the mass parameter for the theory. Instead, 
the mass parameter is determined by t^^'-i^fc^o(*^ ~ ~ ~2/i|. > 0. A characteristic 
feature of SSB is the development of an imaginary sector in Ufj^ki^) in the regime where 
k and $ are small and the argument /c^ -|- p,^ + Xr^^/2 inside the logarithm and the 
square root becomes negative. Any attempt to extract from the low temperature range 
the critical temperature Tc beyond which symmetry restoration takes place using (5.1) is 
fruitless since the resulting Tc is complex due to the presence of an imaginary part in the 
finite temperature contribution found in the last line of (5.1) [4]. However, concentrating 
only on the real part of C//3^/c($) and taking into account the higher order daisy and super- 
daisy graphs, Dolan and Jackiw [4] obtained a critical temperature Tc 



^ ~Xr ^ 



1/2 



(5.3) 



beyond which symmetry restoration takes place via a second order phase transition. Never- 
theless, a true symmetry restoration should be accompanied by the disappearance and not 
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negligence of the imaginary contribution in the thermal blocked potential Up^ki^)- This 
is clearly not possible within the context of the IMA since the finite temperature effects 
never enter the argument /c^ + jl\ + \r^^/2, and we see that the imaginary contribution 
persists for all T unless higher loop effects daisy and super-daisy are properly taken into 
account [4]. 

It is interesting to compare our finite temperature RG prescription with that used 
in [4] since in our approach the argument inside the logarithm and square root becomes 
/c^ + and is a positive quantity beyond T^. 

With the RG flow equation (2.9), we define Tg to be the symmetry restoration tem- 
perature above which the imaginary contribution vanishes entirely. A comparison between 
Tc and is illustrated in Fig. 9. Notice that in the RG scheme, fi'j^ is used instead of 
jlj^ which together with determines Tc- For sufficiently low one finds agreement 

between the RG result and (5.3). However, they begin to deviate for large — /x^ and the 
RG approach gives a higher symmetry restoration temperature compared to Tc in this 
limit. This observation can be explained by noting that for a given T, the RG prescription 
yields a lower value for the thermal mass parameter compared to that of [4] for large T 
due to the fact that we also take into consideration the evolution of the coupling constant 
A^. Therefore, for the condition //^^ = to be satisified, we must have Tc > Tc- Had we 
replaced A/3 by A^ by neglecting the influence of the running of the coupling constant on 
the mass gap equation (4.11), we would have Tc — Tc instead. 

Notice that it is also possible to extract Tc without confronting the complication of 
imaginary contribution by starting from a temperature value T > Tc and gradually lower 
the temperature. Both methods lead to the same T^ as has been checked numerically. 
Near the critical point T > T^ one obtains from (4.11): 

ttT r 1 1 

such that fip^ = at exactly T = Tc- Similarly, a resummation of higher physical loop 
contributions in (4.12) shows that in the vicinity of critical point, the coupling constant 
behaves as: ^ 

Xf3 ^ ^ \l - -4-^ hm {XpT^)] , (5.5) 

after substituting (5.4). One therefore sees that the effective thermal coupling constant Xp 
vanishes at the transition temperature and the theory becomes non-interacting. This can 
also be obtained directly from (4.4) which implies 



31 
8^ 



XI+Xp^-Xr = (5.6) 



at precisely T = T^ where 



Jo z{e--l) ^ ' 

That X diverges naturally gives Xp^ = 0. The temperature dependences of < $ >, and 
A/3 above and below Tc are depicted in Fig. 10. Our numerical result indeed yields A/3^ = 
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when approaching both from above and below. However, for T < in the SSB phase, 
there exists certain numerical fluctuations which nevertheless diminish at Tc. A second 
order phase transition characterized by a continuous decrease of the vacuum condensate 
is predicted with our RG flow equation. With = at the transition point, the IR 
divergence is completely lifted. Had we not incorporated the higher loop effects, the IR 
singularity would persist and result in a breakdown of the perturbation theory. 

The inapplicability of the IMA scheme near Tc can be illustrated by studying the order 
of transition it predicts. For T > Tc, with > and Xp < 0, a flrst order transition is ob- 
tained due to a positive higher order contribution, in contradiction with that predicted 
with RG. Our analyses are in agreements with that obtained in [10] and [17]. However, in 
[10] where a smooth momentum regulator is used, there exists residual dependence of the 
running parameters on the shape of the momentum regularizing function. 

Physically the critical temperature corresponds to a flxed point in the RG trajectory. 
Since this flxed point is of Gaussian nature with fi'^ = Xp^ = 0, any interaction between 
the scalar flelds must be of higher order and parameterized by the irrelevant operators 
classified around the fixed point. For this theory, the critical exponents can be accurately 
determined and shown to coincide with that of the three-dimensional theory at T = [10]. 
From Fig. 10, we also observe that after symmetry restoration with T > Tc, Xp rises again 
and eventually approaches the same constant value as that without going through phase 
transition. 

VI. SUMMARY AND DISCUSSIONS 

In this paper we have carried out the investigation of finite temperature scalar theory 
using an improved RG program. Our RG flow equation has successfully reproduced the 
characteristic behaviors of the system in both high and low T limits. In particular, we have 
established a connection between the 0(4) and 0(3) blocked potentials in large k and small 
k regimes in sec III, and in sec IV dimensional reduction and the proper high temperature 
RG flow patterns for and ^ are deduced. At /c = the coupling constant Xp^k is 
seen to decrease with T and approximately approaches a positive constant value. Thermal 
effects are suppressed for large k and vanish entirely at A; = A where the bare theory is 
defined. In addition we demonstrated in Sec. V how true symmetry restoration is attained 
with vanishing imaginary contribution in Up^ki^)- 

Our method is more advantageous compared to the approaches mentioned in the 
Introduction in a number of ways: The coupled equations (4.9) and (4.10) derived from 
(2.9) for n'p ^ and Xp^k have provided further improvement compared to that derived for 
n'p alone in [4] and [17] since the continuous feedbacks between ^ and A^^^ are retained 
systematically. In addition, through (2.9), we are able to analyze simultaneously the 
competing effects of two scales, namely, k and T in an unambiguous manner. A RG 
approach based on the running of an internal scale k is more "physical" than that employed 
in [6] with the external T as the running parameter. When applied to the Yang-Mills theory 
with a momentum-dependent gluon polarization tensor, the resummation using the latter 
scheme is inconsistent [6]. This poses no difficulty with RG formulation which naturally 
includes both the temperature and momentum effects. 
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Various interesting issues can be explored with our RG scheme in hght of its success. 
For example, one can can use this approach to study the nature of the phase transition 
for the electroweak theory. For A^; << g]^, where qr is the coupling constant for the 
gauge fields, one would expect a first order transition which is required for explaining 
the asymmetry of the baryogenesis [20]. For the Yang- Mills theory, a RG flow equation 
similar to (2.9) will provide information on the roles of k and T on the running of the 
gauge coupling constant. The flow of the theory with at T = has been worked out [21]. 
It would be interesting to investigate the effect of T on such a theory which is known 
to exhibit asymptotic freedom at T = 0. If T and k can generate opposite effects as for 
the scalar theory, there will be nontrivial consequences on the picture of deconfinement 
transition of quarks and gluons. In addition, the resummation of "hot thermal" loops 
using this RG approach will readily yield the gauge-independent gluon damping rate and 
be compared with that obtained in [5] via an effective action. Works along these directions 
are currently in progress. 
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APPENDIX 



We proceed to evaluate the integrals shown in (4.3) and (4.4) following the techniques 
utilized in [4]. The presence of the blocking scale k will slightly complicate the algebra 
compared with the previous studies. We first consider the small a limit of the integral in 
(4.4) having the form 



/(a, 6) 



dz 



^(z2 + a2)(z2 + ^,2) 



(e^^-l) ' 



dz 



^(^2+^2^(^2+52) 



coth 



^ \/2:2 + 



fh: 



2-'0 ^(^2 + ^2) (^2 + ^2) 

= /(i)(a,6) + /(2)(a,6) , 



2V^2 + q2 J2 



1 



z"^ -\- a? -\- 47r2n2 



(A.l) 



where 



^0 



dz 



E 



^2 + ^,2) z^+a^ + 47r2n2 ' 



(A.2) 
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and 



1 



z'^ + a?) (^2 + 



(A.3) 



Since divergences are encountered when splitting the integral in such a manner, we intro- 
duce a suppression factor ^""^ to regularize the individual sum and expect the infinities to 
cancel, thereby making the final result for (A.l) finite. By rewriting (A. 2) as 



OO Icy /.c 

»i— _, — , •'0 



dx 



X 



(1 + x'^)\/h'^x'^ + + Att'^h^ '' 



(A.4) 



via a change of variable x = {a? + A'K'^in?'Y/'^z ^ , the integration can be carried out with 
the help of 



^0 



dx 



X 



(l + a;2)V62x2 + c2 2^ 2 

TT sec(e7r/2) 



+ 



2Vc2 - 62 



where 



F(a, 6, c; 7) = a, c; 7) = 5-^(6, c-h) ! dx x^-^{l - xY'^'^il - 7a;)" 



with 



r(a; + y) Jo 
Substitution of (A. 5) into (A.4) then gives 



(a^ + 47r2n2)' 



(A.5) 
(A.6) 

(A.7) 



/W(a,6) = -sec(-) [ 1^ , T^^'^^^^' 



,1-e 



+ i^r(i-5)r( ^ 



1 + e 



) a-2F(l,l 



e 3 — e 6 



+ 2 J](a2 + 47^V)~'F(l,l- 



2' 2 'a2 

e 3-e 



(A.8) 



which in the limit of vanishing e becomes 

TT 1 



/W(a,6) 



1 ^ 1 



2 - 62 
1 



+ 2-^-%-C(l + e) + -5:^ 

n=l 



2' 2 'a2 + 47r2n2- 

a2_?,2 1/2 



1/2 



- 1 
b 



Va^ - 6^ 

1 TT 

= :^ + 



sm 



n=l * 



1 



2e 2 Va2 - ^2 2 



(A.9) 
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where we have used 



c(i+e) = E 



^'""'^'^^ -(27r)^[^-ln27r + 7 + 0(e)] 



r(l + e) sin(f ) 



oo 

n 

n=l 



2„2 



a + 4-n n 



1 sm ( 



n=l 

= 6 — coth (-) - — r 
L4o ^2^ 2a2 



1 - 



62 n-1/2 



6 5^ (a^ + 47r V) " + 



n=l 



_ b_ 

~ 24 ^ 1440 ^ 



and neglected 



n=l 



2 _ ^2 1/2 



(a' - (>") y. J_ 



n=l 



0{a'). 



In a similar manner, we have 



dz 



^(^2 +q2)(^2 +^2) 



40F\ 



r(i-i)r( 

1 - e,„,e, „A e 1 + e fe^ 



-l + e^ ^/l 



e 3-e ^2 



^^^'^'■^ 2' 2 'a2^ 



2 ' ^2' 

2 

2 '2' ' ^2' 2' 2 ' 



tanh i( )- + ln- + 0(e). 



2 ^a^ 2e 2 2 

Combining (A. 9) and (A. 14) yields the finite result: 



I{a, h) — I dz — = 

•^0 ^(^2 + ^2) (^2 + 52) 

1 

" Va2 - IP 



(e^/^-l)"' 



"TT . -1 




1 r 


— — sm 


(^)] 




L2 







0' , -1 / &\ 

In — + 7 + tanh M-) 

471 ^a' 



b 
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which for vanishing 6, reduces to that obtained in [4] . 
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On the other hand, ior k » T » fj,ji, the original integral expression in (4.4) can 
be approximated as: 



i{a,b) 



dx 



^x'^ - (a2 - 62) - 1 + 



xe 



rj.2 _ 2)2 



1-7(1 + ^)^-:^ 



where we have used 



oo g-a; 



(-1) 



n+1 



n-1 



X 



n+1 



nl 



(-l)'"a'^ 



with 



Ei{-a) = - / — e"^. 



a" n(n — 1) ■ ■ ■ (n — m) 

m=o \ ^ \ I 



dx 



X 



Since Ei{—a) ^ — e "/a as a — >• cxd, we arrive at 



I{a oo, 6) 



1 + a- -(a^ - 6^) + 



The other integral shown in (4.3) has the following form: 



J(a,6) 



Before evaluating (A. 20) fully, we first take the limit 6 = and write 



J(a) = J(a, 0) = / dz 



Following the procedures outlined above, eq. (A. 21) becomes 



2 /» 



z'^ + a'^ + Att"^!!? z^ + 47r2n2 



1 P 
2/0 



(A.16) 



(A.17) 



(A.18) 



(A.19) 



(A.20) 



yz^+a? 



— z 



= ^ + jW(a) + j(^)(a), 
where we have added and subtracted 



roo 00 „ 



+ 47r2 



dz z 



(A.22) 



(A.23) 



TT 
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Using 



^0 



v2-e 



dz 



TT /^^n/ 1-€ l-e\ 



the first integral in (A. 22) in its regularized form can be written as 



jP{a) = -a^ Y: / 



+ 47r2n2)(2;2 + 47r2n2 + a^) 



-f-<=(f) E |(4-V + a^) 



,2^2 , „2\(l-e)/2 



(27rn) 



l-e 



= -f sec(f){a^- + 2 jW(a)} 



where 



oo 



Since 



= a(l - 6) £(47rV + a^)-(^+^^/^ 



n=l 



= a(l-e)< 



C(i + e) 

(27r) 



l+e 



(27rn) 



l+£ 



1 + 



-(l+^)/2 



a(l-6)(27r)-(i+^)C(l + e) + 0(a3), 



this implies 



^i'n«) = y(l-6)(27r)-(^+^)C(l + 6) + 



or, in the vanishing e limit, 

Ji^)(a) = 
In a similar manner, we have 



Tra a 

T ~ T 



- - ln27r + 7 - ij + 



J(^)(a) = -i/ dz 



1 



1 , a 



- -In- 



2J' 



with the help of 

L 
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and discarding the a-independent term. Adding up (A. 29) and (A. 30) then leads the final 
result 



J(a) 



Vz^+a^ 



TT^ 7ra a? 



(A.32) 



a 



6 2 4 
Finally, to evaluate (A.21), we observe that 



In— + 7-1 
47r 



+ 



dJ{a, h) 
db 



bl{a,b). 



(A.33) 



J(a, 6) = !L _ ^ ^/^?^ - 7 (a' - b^) \^^^ + 7 + tanh'^ - 
62 4^ ^L47r 



Integrating over 6 with the help of (A. 15) and imposing the boundary condition (A.32) 
gives 

I/O , o\ r, a ^ -ib~ 

a- 

— ■ -1 ^ , b^ 

a 4 36 

The above expression can be checked by taking the limit a = b where (A.21) can be 
simplified using [22] : 



(A.34) 



du 

Jo e'^-l 



and leads to 



r 



E 2(£+l) ^^(2n + £)(2n) 

(a) = / c 
Jo 



(A.35) 



u 



TT 



1 _ ^ _L V ^^nQ^'" 1 

4^(2n+l)! 

n=l ^ ' 



- 1 



(A.36) 



With i?2 = 1/6, ones finds an amazing agreement between J* (a) and J {a, b — a) in (A.34). 
Similarly, for 6 — > a — > cxo, we have 



J{a,b^ oo) = e " 



a &2 
^+2 + 2^ + 



(A.37) 
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FIGURE CAPTIONS 

Fig. 1. Diagrammatic representation of loop resummation. Contributions from all non- 
overlapping graphs are included in our physical loop. 

Fig. 2. Flow pattern of the blocked potential for T = with = 10~^, = 0.1 and 
A = 10. Solid and dashed lines represent the RG and IMA results, respectively. 
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Fig. 3. Flow pattern of the blocked potential for T = 5 with /ifj = 10 ^, Xr = 0.1 and 
A = 10. Solid and dashed lines represent the RG and IMA results, respectively. 

Fig. 4. Evolution of /i^ ^ and Xp^k as a function of k for various values of T using 
fll^ = 10-4, Xr = 0.1 and A = 20. 

Fig. 5. Comparison of the flow of the mass parameter at T = for different RG prescrip- 
tions using flj^ = 10-4, Xr = 0.1, A = 30, and A = A/y/2. 

Fig. 6. Comparison of the flow of the coupling constant at T = for different RG 
prescriptions using fij^ = 10-^, Xr = 0.1, A = 30, and A = A/-\/2. 

Fig. 7. Limitation on dimensional reduction at high temperature. Notice the gap be- 
tween the mass parameters generated from the dimensionally reduced and the full RG 
prescriptions. 

Fig. 8. Temperature dependence for /j,^ and A/3 at = 0. 

Fig. 9. Comparison of the critical temperature obtained by RG with T^,. Notice that a 
higher value is predicted for RG. 

Fig. 10. Temperature dependence of various quantities in the symmetry broken phase. 
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